Draft version July 27, 2009 

Preprint typeset using 1^1^^ style cmulatcapj v. 08/13/06 



0^ 
O 
O 



6 

(N 
> 
cn 

(N 
(N 

in 
o 
a^ 
o 



FIRST NEUTRINO POINT-SOURCE RESULTS FROM THE 22-STRING ICECUBE DETECTOR 

IceCube Collaboration: R. Abbasi\ Y. Abdou^, M. Ackermann'', J. Adams^, J. Aguilar\ M. Ahlers"', 
K. Andeen\ J. Auffenberg*^, X. Bai^, M. Baker\ S. W. Barwick*^, R. Bay", J. L. Bazo Alba"', K. Beattie^", 
J. J. Beatty^^'^^, S. Bechet^^, J. K. Becker^^, K.-H. Becker'', M. L. Benabderrahmane'*, J. Berdermann^, 
P. Berghaus\ D. Berley^', E. BERNARDiNr"*, D. Bertrand^^, D. Z. Besson^", M. Bissok^', E. Blaufuss^^, 

D. J. BOERSMA^ C. BOHM^**, J. BOLMONT^, S. BOSER^, O. BOTNER^'^, L. BRADLEY^", J. BrAUn\ D. BrEDEr'', 



T. Castermans^\ D. Chirkin\ B. Christy^"', J. Clem^, S. Cohen^^, D. F. Cowen^"^-^^, M. V. D'Agostino" 
M. Danninger^**, C. T. Day^", C. De Clercq^^, L. DEMIR6RS^^ O. Depaepe^", F. DESCAMPS^ P. Desiati^ 
G. de Vries-Uiterweerd^, T. DeYoung^", J. C. Diaz-Velez\ J. Dreyer^"*, J. P. Dumm\ M. R. Duvoort^'', 
W. R. Edwards"', R. Ehrlich^"', J. Eisch\ R. W. Ellsworth"'', O. Engdegard^'', S. Euler", P. A. Evenson^, 
O. Fadiran^", a. R. Fazely", T. FEUSELS^ K. Filimonov", C. Finley\ M. M. Foerster^", B. D. Fox^°, 
A. Franckowiak^*^ , R. Franke^, T. K. Gaisser'', J. Gallagher^", R. Ganugapati\ L. Gerhardt"'-'', L. Gladstone^ 



A. Goldschmidt^", J. A. Goodman^ 



R. Gozzmr'", D. Grant^" 



T. Griesel^", a. Gross"'", S. Grullon\ 



R. M. GuNASiNGHA^^, M. Gurtner", C. Ha^", a. Hallgren"*, F. Halzen\ K. Han^, K. Hanson\ Y. Hasegawa'", 
J. Heise^'', K. Helbing", p. Herquet^\ S. Hickford*, G. C. Hill\ K. D. Hoffman"'', K. Hoshina\ D. Hubert^^, 
W. Huelsnitz""', J. -P. HuLSs", P. O. Hulth^*, K. Hultqvist^*^, S. Hussain^, R. L. Imlay^^, M. Inaba^^, A. Ishihara^ 
J. Jacobsen^ G. S. Japaridze^'', H. Johansson^**, J. M. Joseph^'', K.-H. Kampert", A. Kappes"^^, T. Karg" 
A. Karle\ J. L. Kelley\ P. Kenny^^, J. Kiryluk^ 



G. Kohnen^\ H. Kolanoski^**, L. Kopke^", M. Kowalski^* 



T. Kuwabara', M. 



Labare , S. Lafebre" 



F. KISLAT^ S. R. KLEINl°'^ S. KLEPSER^ S. KnOPS", 

T. KowARiK^", M. Krasberg\ K. Kuehn", 
K. Laihem", H. Landsman^ R. Lauer^, H. Leich^, D. Lennarz^^, 
A. LuCKE^^ J. Lundberg"\ J. LuNEMANN'*", J. Madsen'", P. Majumdar^, R. Maruyama\ K. Mase'", H. S. Matis^", 

C. p. McParland"', K. Meagher^^, M. Merck\ P. Meszaros^^'^", E. Middell^, N. Milke^*, H. Miyamoto^^ 
A. MoHR^^ T. MoNTARULi"''"', R. Morse\ S. M. Movit^'\ K. Munich", R. NAHNHAUER•^ J. W. NAM^ P. Niessen"^, 

D. R. NYGREN"'■^^ S. OdROWSKF", a. OLIVAS^^ M. OlIVO"-', M. ONO^^ S. PANKNIN2^ S. PATTOn"', 

C. Perez de los Heros"', J. Petrovic^^, A. Piegsa^", D. Pieloth^, A. C. Pohl"'''"', R. Porrata'-', N. Potthoff", 
P. B. Price", M. Prikockis^", G. T. Przybylski"', K. Rawlins^^ P. Redl"'', E. Resconi^\ W. Rhode", 
M. RiBORDY^^, A. Rizzo^", J. P. RoDRiGUEsS P. Roth"'', F. Rothmaier"'", C. Rott", C. Roucelle^\ D. Rutledge^", 

D. Ryckbosch^, H.-G. Sander^", S. Sarkar^, K. Satalecka^, S. Schlenstedt^, T. Schmidt^^, D. Schneider^ 
A. ScHUKRAFT^'', O. ScHULZ^^ M. ScHUNCK^'', D. Seckel', B. Semburg^', S. H. Seo^*^, Y. Sestayo^^ S. Seunarine"*, 

A. SiLVESTRi**, A. Slipak °, G. M. Spiczak^'', C. Spiering^, M. Stamatikos", T. Stanev^, G. Stephens^", 
T. Stezelberger^", R. G. Stokstad^", M. C. Stoufer"', S. Stoyanov^, E. A. Strahler\ T. Straszheim"'', 
K.-H. SuLANKE'^ G. W. SuLLivAN^^ Q. SwiLLENs", I. Taboada''**^ , O. Tarasova^, A. Tepe'*, S. Ter-Antonyan^^, 

C. Terranova^^ S. Tilav'', M. Tluczykont^, P. A. Toale^", D. Tosr\ D. Turcan^^, N. van Eijndhoven^'"', 
J. Vandenbroucke'-', a. Van Overloop^, B. Voigt^, C. Walck^*, T. Waldenmaier^*, M. Walter^, C. Wendt\ 
S. Westerhoff\ N. Whitehorn\ C. H. Wiebusch^', A. Wiedemann", G. Wikstrom^*^, D. R. Williams^'"", 
R. Wischnewski^, H. Wissing^^'^^, K. Woschnagg", X. W. Xu^^, G. Yodh**, and S. Yoshida''^^ 

Draft version July 27, 2009 

ABSTRACT 

We present new results of searches for neutrino point sources in the northern sky, using data recorded 
in 2007-08 with 22 strings of the IceCube detector (approximately one-fourth of the planned total) 
and 275.7 days of livetime. The final sample of 5114 neutrino candidate events agrees well with the 
expected background of atmospheric muon neutrinos and a small component of atmospheric muons. 
No evidence of a point source is found, with the most significant excess of events in the sky at 2.2 cr 
after accounting for all trials. The average upper limit over the northern sky for point sources of 
muon-neutrinos with E""^ spectrum is < 1.4 x 10"^^ TeVcm"'^ s^^, in the energy range from 

3 TeV to 3 PeV, improving the previous best average upper limit by the AMANDA-II detector by a 
factor of two. 

Subject headings: acceleration of particles — cosmic rays — neutrinos 
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1. INTRODUCTION 

Cosmic rays with energies up to lO'^" eV pervade the 
Universe, but their sources remain unknown. Possible 
acceleration sites of cosmic rays include shock fronts 
in supernova remnants, pulsars, microquasars, active 
galactic nuclei, and gamma-ray bursts. While many of 
these sources are n ow observed by g a mma-ray astron- 
omy experiments (lAbdo et al.l l2007t lAharonian et al.l 
120061: [Albert et all [2006D, it remains difficult to deter- 
mine whether the gamma-ray emission is of leptonic or 
hadronic origin. Hadronic acceleration in these sources is 
expected to produce a correlated neutrino flux as accel- 
erated protons interact with ambient gas and radiation 
to produce mesons, and the charged r nesons decay to 
neutri n os (for reviews see iBeckeil ()2008[ ) : iBednarek et al.l 
()2005f ): [Halzen fc Hooper! ( 2002f l and references therein). 
This signature uniquely distinguishes hadronic from lep- 
tonic processes, and thus detection of the high-energy 
neutrino flux from cosmic ray accelerators is the key to 
identifying them. Moreover, since neutrinos can propa- 
gate freely through dense environments and across cos- 
mological distances that are optically thick to photons, 
they can probe hidden regions and reveal unexpected 
sources, opening a unique window on the high-energy 
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Universe. 

Previous searches for high-energy as trophysical neu- 
trino sources performed by MACRO dA mbrosi o et abl 
I2001h . Sup er-Kamiokande (iDesai et al J 2008), and 
AMANDA (jAbbasi et "all 12009a') have set upper Hmits 
that demonstrate the need for much larger experiments 
to detect these weakly interacting particles. The Ice- 
Cube Neutrino Observatory is now instrumenting a cu- 
bic kilometer of the clear Antarctic ice sheet at the geo- 
graphical South Pole. Construction began in the austral 
summer 2004-05, and is planned to finish in 2011. The 
full detector will comprise 4,800 Digital Optical Mod- 
ules (DOMs) deployed on 80 strings between 1.5-2.5 km 
deep within the ice, a surface array (IceTop) for observ- 
ing extensive air showers of cosmic rays, and an addi- 
tional dense subarray (DeepCore) in the detector cen- 
ter for enhanced low-energy sensitivity. Each DOM con- 
sists of a 25 cm diameter Hamamatsu ph oto-multiplier 
tube, e lectronics for waveform digitization (jAbbasi et al.l 
l2009bf ). and a spherical, pressure-resistant glass hous- 
ing. Waveforms are recorded when nearest or next-to- 
nearest DOMs fire within ± 1 microsecond; event trig- 
gers occur when eight DOMs record waveforms within 
5 microseconds. Calibrations ensuring nanosecond tim - 
ing precision are described in lAchterberg et al] (120061 ) . 
The DOMs detect Cherenkov photons emitted by rel- 
ativistic charged particles passing through the ice. In 
particular, the directions of muons (either from cosmic 
ray showers above the surface, or neutrino interactions 
within the ice or bedrock) can be well reconstructed from 
the track-like pattern and timing of hit DOMs. Identifi- 
cation of neutrino-indu ced muon events in Ice Cube has 
been demonstrated in ,Achterberg et aTl ()2007a[ ) using at- 
mospheric neutrinos as a calibration tool. 

2. EVENT SELECTION AND ANALYSIS 

As of spring 2007, there were 22 deployed IceCube 
strings. The physics run for the 22-string configuration 
started 2007 May 31 and ended 2008 April 4, when the 
40-string configuration began operating. The final live- 
time is 275.7 days, about 90% of the total available time 
including operation during the construction season. The 
event trigger rate is ~ 550 Hz, predominantly due to 
down-going muons. The rate of atmospheric neutrino- 
induced muons triggering the detector is roughly 10^ 
times lower. Because only neutrino-induced muons can 
travel upwards, neutrino events can be isolated by se- 
lecting up-going tracks. An online event filter makes 
the first rejection of down-going tracks. Events which 
pass (at a rate ~ 20 Hz) are sent over satellite to the 
North. Likelihood-based track reconstructions are per- 
formed, improving the directional accuracy and back- 
ground rejection capabilities as well as providing ind i- 
vidual angular uncertainty estimates ()Neunhoffeill2006D . 
The reduced log-likelihood of the best-fit track, the angu- 
lar uncertainty, and the number of modules which were 
hit by direct Cherenkov photons (within a window of —15 
to -|-75 ns for estimated hit times of the reconstructed 
track) are the main parameters used to select up-going 
neutrino candidates and reject background. An addi- 
tional cut on the likelihood ratio of the best-fit track to 
the best-fit track constrained to be down-going further 
reduces background close to the horizon. In the final 
analysis, a wide range of cuts based on these parameters 
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are compatible with optimal sensitivity to both hard and 
soft spectrum sources. Within this range, applying cuts 
that also remove the largest fraction of mis-reconstructed 
down-going events yields a final sample consisting of 5114 
neutrino candidate events. 

At this time there is not a single strategy for point 
source searches in IceCube, so more than one approach 
was investiga ted, includ ing a binned analysis similar to 
that in Achter berg et al . ( 2007b ) , as well as an u nbinned 
analysis similar to that in lAbbasi et al.l ()2009aD . Simu- 
lation studies using the 22-string IceCube configuration 
showed that the latter approach was on average 35% 
more sensitive for both hard and soft point source spec- 
tra. For this reason, it was decided before the data were 
unblinded that the unbinned analysis would be used for 
the final results, which are reported below. 

The unbinned like lihood analysis is described in detail 
in iBraun et al.l (pOOSi) : it uses both the direction and 
energy information of each event. Astrophysical neu- 
trino source spectra are typically expected to be harder 
(~ in the Fermi model of cosmic ray acceleration) 
than the known spectrum of the atmospheric neutrino 
background (~ E~^ ''). Thus a neutrino point source 
may be detectable not just by the clustering of event ar- 
rival directions, but by a different event energy distribu- 
tion than the background. For each direction in the sky 
tested, the analysis performs a fit for the number of sig- 
nal events Us above background, and the spectral index 
7 of the excess events. The test-statistic in the analysis 
is the log likelihood ratio of the signal hypothesis with 
best-fit parameters (ng and 7) to the null hypothesis of 
no signal present [us = 0). This test-statistic provides 
an estimate of the significance (pre-trial p- value) of devi- 
ation from background at a given position in the sky. As 
described below, the post-trial significance is determined 
by applying the analysis to scrambled data sets, in which 
the right ascension of the events are randomized but all 
other event properties are kept the same. 

Two unbinned point-source searches are performed. 
The first is an all-sky search within the declination range 
—5° to -1-85°: the maximum likelihood ratio is evaluated 
for each direction in the sky in steps of 0.25° r.a. and 
0.25° dec, (well below the angular resolution of 1.5°). 
The significance of any spot is given by the fraction of 
scrambled data sets containing at least one spot with a 
log likelihood ratio higher than the one observed in the 
real data. This fraction is the post-trial p-value. Because 
the all-sky search involves a large number of effective tri- 
als, the second search is restricted to the directions of 28 a 
priori selected source candidates, in order to improve the 
confidence of a possible detection of one of these objects. 
The post-trial p-value is again found by performing the 
source list analysis on scrambled data sets. The smallest 
post-trial p-value from either of the two searches is then 
taken as the final significance of the analysis, with a final 
trial factor of two. 

3. DETECTOR RESPONSE 

A simulation of and was used to determine the 
effective area and point spread function for the 22-string 
IceCube configuration, shown in Fig. [T] The sky- 
averaged median angular reconstruction error is 1.5° for 
both E^'^ and atmospheric spectra. For an E^"^ (atmo- 
spheric) neutrino spectrum, 90% of the events are in the 




Fig. 1. — Solid-angle-averaged effective areas at final cut level 
for astrophysical muon neutrino fluxes (v^ + at different decli- 
nations. The turnover at high-energy for up-going events is due to 
absorption by Earth. Inset: E~'^ point spread function (angular 
difference between the neutrino and reconstructed muon track) for 
the same declination ranges. 
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Fig. 2. — Sensitivity to a point-source flux as a func- 

tion of decli nation, shown for the final AMANDA-II analysis 
ijAbbasi et al. 2009a), the current twenty-two string IceCube anal- 
ysis, and the predicted sensitivities for the ANTARES experiment 
llAguilar Sancheai2007l) and the final IceCube configuration. 

central range 3TeV - 3PeV (250 GeV - 16TeV). The 
sensitivity (median upper limit , follow ing the ordering 
principle of iFeldman fc Cousins! (|1998D ) as a function of 
declination is shown in Fig. O For a source at declina- 
tion 5° with E^"^ spectrum, approximately 13 (16) signal 
events are needed for a 50% chance of a post-trial 5 tr de- 
tection based on the source list search (all sky search). 

Systematic uncertainties arise chiefly from the model- 
ing of light propagation in the ice. Simulation of muons 
passin g through the ice is described in Chirkin & Rhodi 
(|2004f l. and the propagation of light from the muon 
to t he optical modules is performed with Photon- 
ics (jLundberg et al.l l2007l ). Embedded impurities in 
the ice alter the scattering and absorption properties 
with a depth dependence due to climate changes over 
tens of thousands of years; direct measurements of 
these ice properties are used in the detector simulation 
(jAckermann et akl l2006l ). Data-Monte Carlo compar- 
isons using both trigger level and high-quality down- 
going muons indicate the systematic uncertainty which 
remains from this aspect of simulation is 15%. System- 
atic uncertainties on the optical module efficiency and 
true point spread function contribute an additional 9% 
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Fig. 3. — Declination distribution of the final data sample and 
simulated atmospheric neutrino events with 30% theoretical uncer- 
tainty on th e absolute normalization of the predicted neutrino fiux 
HBari- et al.ll2004 . l200i) . 

uncertainty on neutrino flux estimates, while theoretical 
uncertainties on muon energy losses and th e ncutrino- 
nucleon cross-section (jPumplin et all l2002t l contribute 
5%. The total estimated systematic uncertainty of 
18% has been incorporated into the sensitivity and up- 
per limit calc ulations u sing the m ethods described in 
IConrad et all |2003) andlM] Hool). Because taus from 
tau neutrinos also decay to muons (with a 17% branch- 
ing ratio) , any assumption of an additional component of 
tau neutrinos in the astrophysical flux arriving at earth 
would lead to a tighter upper limit on the muon neutrino 
flux. 

Atmospheric neutrino events, while posing a substan- 
tial background to extra-terrestrial neutrino searches, 
provide a useful verification of the detector and simu- 
lation. In Fig. [3l the declination distribution is shown 
for the data and atmospheric neutrino simulation with 
final cuts applied. The atmospheric neutrino model has 
a theoretical uncertainty of 30 % in the flux normaliza- 
tion in the TeV energy range ()Barr et al.l |2004 120061 : 
lAchterberg et al.ll2007a[) . This uncertainty does not af- 
fect astrophysical flux calculations, but limits the preci- 
sion of atmospheric neutrinos as a check of the detector 
simulation. In 275.7 days of livetime, 4600 ± 1400 at- 
mospheric neutrinos are expected, along with an addi- 
tional component of 400 ± 200 mis-reconstructed down- 
going muons frorti air showers (simulated with CORSIKA 
(|Heck et al.lll998[ )). mainly near the horizon. Within un- 
certainties, this is in agreement with the 5114 events ob- 
served. 

4. RESULTS 

The results of the all-sky search are shown in Fig. [4l 
The most significant deviation from background is lo- 
cated at 153.4° r.a., 11.4° dec. The best-fit parameters 
are = 7.7 signal events above background, with spec- 
tral index 7 = —1.65. The pre-trial estimated p-value 
of the maximum log likelihood ratio at this location is 
7 X 10~^. The post-trial p-value is determined by per- 
forming the analysis with the right ascension of the data 
randomized: 67 out of 10,000 scrambled data sets yielded 
a more significant excess somewhere in the sky. 

The results of the point-source search in the direction 
of 28 source candidates selected a priori are given in Ta- 
ble [TJ The smallest pre-trial estimated p-value is 7% (for 
the TeV blazar lES 1959+650); 66 out of 100 scrambled 
data sets have a more significant excess for at least one 



TABLE 1 

Results for the a priori SOURCE candidate list 



Object 


r.a. [°] 


dec. [°] 


n 


S 


P 


*90 


B20 


MGRO J2019+37 


304.83 


36.83 


3 


1 


0.25 


25.23 


2.6 


MGRO J1908+06 


287.27 


6.28 








- 


7.06 


3.7 


Cyg OB2 


308.08 


41.51 








- 


15.28 


2.6 


SS 433 


287.96 


4.98 


2 


8 


0.32 


11.65 


4.1 


Cyg X-1 


299.59 


35.20 








- 


14.60 


2.4 


LS I -1-61 303 


40.13 


61.23 








- 


22.00 


3.0 


GRS 1915+105 


288.80 


10.95 








- 


7.64 


3.3 


XTE J1118+480 


169.54 


48.04 


2 


5 


0.082 


40.62 


2.8 


GRO J0422+32 


65.43 


32.91 








— 


14.10 


2.2 


Geminga 


98.48 


17.77 










9.67 


2.6 


Crab Nebula 


83.63 


22.01 








- 


10.35 


2.4 


Gas A 


350.85 


58.81 








- 


20.22 


3.5 


Mrk 421 


166.11 


38.21 








- 


14.35 


2.8 


Mrk 501 


253.47 


39.76 








- 


14.44 


2.7 


lES 1959+650 


300.00 


65.15 


5 





0.071 


59.00 


3.2 


lES 2344+514 


356.77 


51.70 








- 


17.94 


2.8 


H 1426+428 


217.14 


42.67 








- 


15.64 


2.7 


lES 0229+200 


38.20 


20.29 










10.24 


2.4 


BL Lac 


330.68 


42.28 


1 


6 


0.37 


22.81 


2.7 


S5 0716+71 


110.47 


71.34 


1 


9 


0.31 


44.76 


3.3 


3G66A 


35.67 


43.03 


2 





0.31 


25.70 


2.8 


3C 454.3 


343.49 


16.15 










9.07 


2.6 


4C 38.41 


248.82 


38.13 










14.29 


2.8 


PKS 0528+134 


82.74 


13.53 










8.27 


3.2 


3C 273 


187.28 


2.05 





9 


0.37 


11.73 


4.1 


M87 


187.71 


12.39 










7.91 


3.2 


NGC 1275 


49.95 


41.51 


2 


2 


0.21 


28.32 


2.6 


Cyg A 


299.87 


40.73 










15.05 


2.6 



Note. — Us is the best-fit number of signal events; when rts > 
the (pre-trial) p-value is also calculated. 'I'go is the upper limit of 
the Feldman-Cousins 90% confidence interval for an fiux, i.e.: 
d'S>/dE < $90 10-i2rpgY-lpjjj-2g-i(£;yrpgyj_2 -pj^g backgrouud 

event density at the source declination is indicated by the mean 
number of background events B20 expected in a bin of radius 2° . 

source on this list. 

Of the two searches, the most significant result comes 
from the all-sky search. Accounting for this last trial 
factor of two, the final p-value for the analysis is 1.34%. 
At this level of significance, the excess is consistent with 
the background-only null hypothesis. If not a statisti- 
cal fluctuation, the excess will be detectable with future 
IceCube data, unless it were caused by a one-time or 
rare astronomical event. Subsequent examination of the 
times of the events in the region of excess, however, has 
not revealed any burst-like distribution in time, with the 
ten events that contribute most to the excess distributed 
throughout the year and each separated by a minimum 
of nine days from the next. 

5. GONGLUSIONS 

A search for point sources of high-energy neutrinos has 
been performed using data recorded during 2007-08 with 
22 strings of IceCube. An all-sky search within the dec- 
lination range —5° to +85° found the most signiflcant 
deviation from the background at 153.4° r.a., 11.4° dec. 
Accounting for all trials in the point source search, the 
flnal p-value for this result is 1.34%, consistent with the 
null hypothesis of background-only events at the 2.2 a 
level (if the p-value is expressed as the one-sided tail of 
a Gaussian distribution). No obvious source candidates 
are near this location, and an analysis of the timing of the 
events did not flnd any evidence of a burst in time. The 
location can be added to the a priori source candidate 
list for analysis using future IceCube data, in which case 
a similar excess would be identified with much higher 
significance. 
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Fig. 4. — Equatorial skymap of events (points) and pre-trial significances (p-value) of the all-sky point source search. The solid curve is 
the galactic plane. 



No evidence of neutrino emission was found for any 
of the 28 a priori selected point-source candidates, and 
the resulting upper limits severely constrain some mod- 
els of neutrino emission. For example, simulation of 
the neutrino flux from the r nicroq uasar SS433, using 
the model of iDistefano et al.l (12002 ) and expressed as a 
broken power-law ("Be dnarek et aL [2005) with a cutoff 
at 100 TeV, predicts a point-source signal of 48 neu- 
trino events in the 22-string data sample, and is now 
well excluded. Alte rnative predictions for SS433 (e.g. 
iRevnoso et al.l [20081 ) will however only be testable with 
the full IceCube detector. Correlation of the neutrino 
arrival times with the known orbital period of the binary 
system or the precession period of the jets may enhance 
the sensitivity to such objects and provide insight into 
the emission processes. 

The sensitivity of this search with one season of 22- 
string data already exceeds the combined sensitivity of 
all previous neutrino point-source searches in the TeV- 
PeV energy range. New searches are underway to ex- 
tend the sensitivity to ultrahigh-energy sources in the 
southern sky, and to lower energy sources using events 
recorded by the combined IceCube-AMANDA detector. 
With completion of the full 80-string detector expected 
in 2011, the improved acceptance, signal efficiency, back- 
ground rejection and angular resolution (< 0.8°) should 



provide more than an order of magnitude enhancement 
in sensitivity within several years of operation. 



We acknowledge the support from the following agen- 
cies: U.S. National Science Foundation-Ofhce of Polar 
Program, U.S. National Science Foundation-Physics Di- 
vision, University of Wisconsin Alumni Research Foun- 
dation, U.S. Department of Energy, and National Energy 
Research Scientific Computing Center, the Louisiana 
Optical Network Initiative (LONI) grid computing re- 
sources; Swedish Research Council, Swedish Polar Re- 
search Secretariat, and Knut and Alice Wallenberg Foun- 
dation, Sweden; German Ministry for Education and 
Research (BMBF), Deutsche Forschungsgemeinschaft 
(DFG), Germany; Fund for Scientific Research (FNRS- 
FWO), Flanders Institute to encourage scientific and 
technological research in industry (IWT), Belgian Fed- 
eral Science Policy Office (Belspo); the Netherlands Or- 
ganisation for Scientific Research (NWO); M. Ribordy 
acknowledges the support of the SNF (Switzerland); 
A. Kappes and A. Grofi acknowledge support by the EU 
Marie Curie OIF Program; J. P. Rodrigues acknowledge 
support by the Capes Foundation, Ministry of Education 
of Brazil. 



REFERENCES 



Abbasi, R., et al. 2009a, Phys. Rev., D79, 062001 
— . 2009b, Nucl. lustrum. Meth., A601, 294 
Abdo, A. A., et al. 2007, Astrophys. J., 664, L91 
Achterberg, A., et al. 2006, Astropart. Phys., 26, 155 
— . 2007a, Phys. Rev., D76, 027101 
— . 2007b, Phys. Rev., D75, 102001 

Ackermann, M., et al. 2006, J. Geophys. Res., Ill, D13203 

Aguilar Sanchez, J. A. 2007, in Proc. 30th ICRC, Merida 

Aharonian, F., et al. 2006, Astrophys. J., 636, 777 

Albert, J., et al. 2006, Science, 312, 1771 

Ambrosio, M., et al. 2001, Astrophys. J., 546, 1038 

Barr, G. D., Gaisser, T. K., Robbins, S., & Stanev, T. 2006, Phys. 

Rev., D74, 094009 
Barr, G. D., et al. 2004, Phys. Rev., D70, 023006 
Becker, J. K. 2008, Phys. Rept., 458, 173 

Bednarek, W., Burgio, G. F., & Montaruli, T. 2005, New Astron. 
Rev., 49, 1 

Braun, J., et al. 2008, Astropart. Phys., 29, 299 



Chirkin, D., & Rhode, W. 2004, hep-ph/0407075 

Conrad, J., Botner, O., Hallgren, A., & Perez de los Heros, C. 2003, 

Phys. Rev., D67, 012002 
Desai, S., et al. 2008, Astropart. Phys., 29, 42 

Distefano, C, Guetta, D., Waxman, E., & Levinson, A. 2002, 

Astrophys. J., 575, 378 
Feldman, G. J., & Cousins, R. D. 1998, Phys. Rev., D57, 3873 
Halzen, F., & Hooper, D. 2002, Rept. Prog. Phys., 65, 1025 
Heck, D., et al. 1998, FZKA-6019 
Hill, G. C. 2003, Phys. Rev., D67, 118101 
Lundberg, J., et al. 2007, Nucl. lustrum. Meth., A581, 619 
Neunhoffer, T. 2006, Astropart. Phys., 25, 220 
Pumplin, J., et al. 2002, JHEP, 07, 012 

Reynoso, M. M., Romero, G. E., & Christiansen, H. R. 2008, Mon. 
Not. Roy Astron. Soc, 387, 1745 



